Electron scattering due to dislocation wall strain field in GaN layers 
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The effect of edge-type dislocation wall strain field on the Hall mobility in n-type epitaxial GaN 
was theoretically investigated through the deformation potential within the relaxation time ap- 
proximation. It was found that this channel of scattering can play a considerable role in the low- 
temperature transport at the certain set of the model parameters. The low temperature experimental 
data were fitted by including this mechanism of scattering along with ionized impurity and charge 
dislocation ones. 

I. INTRODUCTION 

As is known GaN films are under extensive examination for many years because of their promising application for 
the construction of short-wavelength light emitting devices^. However, the performance of these devices is limited by 
defects, both native and impurity types. Native defects, in particular, are threading dislocations with high densities 
(10 8 — 10 11 cm~ 2 ) which are result from the large lattice mismatch between epilayer and substrate^! 5 -. Dislocations, 
being charged objects, act as scattering centers (core effect) for carriers affecting the transverse mobility in films 6 . 
In most studies in the context of the GaN layers this Coulomb scattering has been consideredi 8 -^^. At the same 
time, in addition to the core scattering, dislocations can give contribution to the resistivity through deformation and 
piezoelectric potentials&ii. In GaN layers these potentials associate with the strain field of dislocation arrays which 
form low-angle grain boundaries or separate dislocations in the specific casesi. However, the carrier scattering due to 
piezoelectric potential has been found as negligibly small within the bulk of GaNi^. 

II. MODEL 

In this paper we theoretically investigate the contribution to the Hall mobility in GaN layers from a wall of 
dislocations of edge type. The scattering of electrons by dislocation wall (DW) is treated in the framework of the 
deformation potential approach. In this case, the perturbation energy of electron can be written in the standard 
form 1 ^ 

SU(r) = GA(r), (1) 

where A(r) is the dilatation of lattice around a finite dislocation array, G is the deformation-potential constant. 

Let the threading dislocation segments with coordinates (0, h) along the (0, 0, 1) axis (perpendicular to the in- 
terface/layer plane) form a dislocation wall of finite length 2L. The dilatation at the point r > 2L around such 
defect, then, can be found based on the so-called disclination model of grain boundaries and dislocations for isotropic 
mediu m 15 i 16 . It takes the formii 
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where p\ = (x ± L) 2 + y 2 , b is the Burgers' vector of dislocation in the wall, p is the distance between dislocations 
in the wall, a is the Poisson constant. It is interesting to note that Eq.(2) is also the exact formula for the dilatation 
around the high-angle grain boundary as well as the disclination dipole. 

The square of the matrix element of electron scattering in momentum states k to states k with the perturbation 
energy given by Eq.(2) deduced ir*i£ as 
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where A — ^ — \. As one can see from Eq.(2) and Eq.(3), the scattering due to perturbation SU(r) is three 
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dimensional, therefore (k — k ) z = q z = <J 'q 2 — q\ ^ in contrast to the case of an infinitely long dislocation line 
along the (0,0, l)-axis. Omitting the details of further calculations, we can come to the following equation for the 
relaxation time due to the strain field of DW 
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(j 2 (2Lfc ± ) + J 2 (2Lk ± ) - 2^-J (2L*±)Ji(2LA!±)), (4) 

where n^ef is the arial density of DW which is inversely proportional to the square of the mean distance between 
grain boundaries, m* is the effective mass of carrier, J n (t) are the Bessel functions, Si(x) is the sine integral function, 
k±_ = (k x , k y ) is the normal to the disclination line component of the wave vector, h is the reduced Planck's constant. 
When a film is thick (k±h >> 1, bulk regime) Eq.(4) can be simplified to 
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(j 2 (2Lfc ± ) + J*{2Lk ± ) - -L-J^Lk^J^Lk^, (5) 
The relaxation time due to Coulomb scattering at charged dislocation lines can be written in the well-known form 
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T*is(k) = - ^ , (6) 



where N^is is the dislocation density in the wall, c is the distance between acceptor centers along the dislocation line, 

/ is the occupation rate of the acceptor centers along the dislocation, e is the dielectric constant, e is the electronic 
fck T 

charge, \d = y — % — 1S * ne Debye screening length with electron concentration n = — iVJ — f{Ndi S /c), and ks 
is the Boltzmann's constant. To determine the filling factor /, the procedure fromiS has been used. 

Assuming the validity of the non-degenerate statistics, we can evaluate the DW contribution /i wa u to the total 
mobility using the well-known formula 

(t\ - ^ 2 Jr(k)k 2 f d 3 k _ e(r) 
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where fcj = k x r y -\ is the planar component of the wave vector, and fo is the Boltzmann distribution function. 

III. NUMERICAL RESULTS 

We first analyzied the drift mobility contribution due to grain boundary strain field together with other important 
mechanisms of scattering at low temperatures (ionized impurities, charged dislocation lines). The component of the 
scattering due to ionized impurities has been taken in the form of the Brooks-Herring-like formula obtained on the 
basis of the partial- wave phase-shift method 19 . This approach yields the correct results for low T and high n where 
the Born approximation can be false. It takes the form 

-^— = N I va B (k)H , (8) 

where N] is the ionized impurity concentration, v is the electron velocity, cr B (k) is the Born cross section, and Hq is the 
factor of correction obtained from the phase-shift calculations^ 9 -. We found that the Coulomb dislocation scattering 
is major at free carriers concentration n < 10 17 cm~ 3 when Ndis ~ 10 8 cm~ 2 (T < 100K), and above this n the 
impurity scattering dominates. In our calculations dislocation core scattering always dominates above Ndi S ~ 10 s 
cm~ 2 . The numerically calculated fi wa ii on the basis of the Eqs.(4),(6) is shown in Fig.l as a function of temperature 
for some selected model parameters together with other contributions. The deformation-potential constant G has 
been taken equal to 4 eV, that corresponds to the typical values for semiconductors. From our analysis we found 
that ix W aii ~ T 3 / 2 (unlike the case of the separate dislocations where fj, ~ T), and this contribution can be maximal 
for some chosen set of the parameters when the concentration of the DW is sufficiently high (ridef — 10 10 cm -2 ) (see 
Fig.l). 

The low temperature part of the experimental data from^I for two GaN samples with Ndi S = 4 x 10 s cm -2 and 
2 x 10 10 cm~ 2 on the sapphire substrate has been fitted based on the formula for the Hall mobility fin = e(r 2 )/m*(r). 
Here the total relaxation time is given by 

wall ' (lis ' ii 

The results are presented in Fig. 2. As shown, there is agreement with the experimental data when 30K< T < 100K 
for both samples. We found that the result of the fit essentially depends on the distance between dislocations in the 
wall p, as compared to other parameters, and, hence on the angle of misorientation between grains (9 = arcsin(fe/2p)). 

Our preliminary results show that this effect of the angle variation on the Hall mobility can be noticeable even for 
a narrow interval of 9 between 1° and 5°. The second very sensitive parameter in our calculations is the position of 
the dislocation acceptor level energy referred to the conduction band edge energy. 

As noted above, Eq.(2) describes both the dilatation around the low-angle grain boundary and disclination dipole. 
The concept of disclination dipole has been applied to obtain the deformations around high-angle grain boundaries 
and linear defects of the rotational type^. In this connection, the observed in hexagonal GaN layers deformations 
associated with the 5/7 and 4/8 rings and high-angle grain boundaries (see, for example, Refsj 23 ) 24 ) can be considered 
in the framework of the disclination dipole model. The second point which should be noted concerns the interfacial 
misfit dislocations 25 . Misfit dislocation scattering along with considered in this paper can be given in the framework 
of the multi-layer model proposed in 2 ^. 

IV. CONCLUSIONS 

In this article we have theoretically investigated the possible role of the strain field associated with the dislocation 
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wall on the mobility in hexagonal GaN layers. It has been found that this contribution to the total transverse mobility 
can be noticeable at low temperatures for given above densities of such defects and deformation constant typical for 
semiconductors. Our calculations show the core scattering due to charged dislocation lines is dominant mechanism 
when Ndis > 10 9 cm~ 2 . This supports findings in previous publications devoted to this material. At lower dislocation 
densities the ionized impurity scattering and strain field scattering can dominate. 
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Figure Captions 



Fig.l. Calculated contributions to the total drift mobility as a function of temperature for the model parameters: 
N dls = 10 8 cm~ 2 , n def = 1.5 x 10 10 cm" 2 , / = 0.8, n=4 x 10 16 cm" 2 , p = 10" 2 /im, h = 1.2/im. The rest of the 
parameter set has been taken from Refj^i 



Fig. 2. Hall mobility vs temperature for two samples with Ndi S — 4 x 10 8 cm 2 (squares), N d is = 2 x 10 10 cm 2 (circles) 
from Ref. Solids lines are theoretical curves. The set of the model parameters: (for squares) G — 9 eV, / = 0.92, 
p = 9x 10"Vm, n def = 2.7 x 10 10 cm" 2 , (for circles) G = 7 eV, / = 0.98, p = 6x 10"Vm, n def = 1.5 x 10 10 cm" 2 . 
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